Investigation of the seismicity after the initiation of a Seismic Electric Signal activity 

until the main shock 
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The behavior of seismicity in the area candidate to suffer a main shock is investigated after the 
observation of the Seismic Electric Signal activity until the impending mainshock. This makes use 
of the concept of natural time \ an d reveals that the probability density function of the variance 
k\{— {x 2 } — (x) 2 ) exhibits distinct features before the occurrence of the mainshock. Examples are 
presented, which refer to magnitude class 6.0 earthquakes that occurred in Greece during the first 
two months in 2008. 
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I. INTRODUCTION 

Seismic Electric Signals (SES) are transient low fre- 
quency (< 1Hz) electric signals that have been observed 

in GreecedHlalllEIEIa, Ja P an B[l, Mexico^ etc. 
days to months before earthquakes (EQs). They are emit- 
ted when the stress in the focal region reaches a critical 
value before the failure [HI, G3- This stems from the fact 
that a stress variation affects the Gibbs energy for the 
defect formation [l3j], migration[l4| and activation[l5| in 
solids. For EQs with magnitude 6.5 or larger, SES are ac- 
companied by detectable magnetic field variations [l6l Il7l 
[l8j . A sequence of SES observed within a short time (e.g. 
« lh) is termed SES activity the analysis of which has 
been shown to obey an 1/ f -behavior [l9l [20| . Recently, 
a method [2l|, |22J, |23j, \2A\, |25j, |26| has been presented that 
enables the shortening of the time- window of the impend- 
ing mainshock from a few hours to a few days only. It is 
based on the concept of a new time domain, termed natu- 
ral time E3, El 



22| and investigates the order parameter 



of seismicity [23j (see also below) that occurs after the 
SES activity and before the main shock in the area can- 
didate to suffer a strong EQ. The improvement of that 
method constitutes the basic aim of the present paper in 
view of the great practical importance in determining the 
time of an impending catastrophe. Along these lines, the 
most recent SES electric field data are also presented. 

In a time series consisting of N events, the natural 
time Xk = k/N serves as an index (TqI. El E3 for the 
occurrence of the fc-th event. The evolution of the pair 

hk^Qkijs studied [H, m, El El, El El, ES EE El El, 

|29|. |30L l3l|. where Qk denotes a quantity proportional to 
the energy released in the fc-th event. For dichotomous 
signals, for example, which is frequently the case of SES 
activities, Qk can be replaced by the duration of the fc-th 
pulse. As a second example, we refer to the analysis of 
seismicity [H El El, El], where Qk may be considered 
as the seismic moment Mqu of the fc-th event, since Mq 
is roughly proportional to the energy released during an 



EQ. The normalized power spectrum is given [1 
by II(w) = |<&(a;)| 2 , where $(w) is defined as 



21, 22 
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p fe exp | iuj— 



(1) 



In this definition, p k stands for p k = Q k /J2n=iQn, 
and lj = 27T0; where </> denotes the natural frequency. 
The continuous function in Eq.{T]) should not be 

confused with the usual discrete Fourier transform be- 
cause the latter considers only the relevant values at 
4> = 0, 1, 2, . . ., while in natural time analysis the proper- 
ties of II(w) or II(</>) are studied[H, [lj, El El for nat- 
ural frequencies <f> less than 0.5. This is so, because in 
this range of <j), H(u>) or H((f)) reduces to a characteristic 
function for the probability distribution p k in the context 
of probability theory. 

When the system enters the critical stage, the following 
relation holds [H El El : 



6cosa-> 12 sin w 



5lu 2 



5w 3 



(2) 



which for w — ► 0, simplifies to|13. Il9l. l2l| 

IL(u) w 1 - 0.07w 2 . 

This relation reflects 23] that the variance (x 2 ) ~ (x) 2 °f 
X is given by 



(X 2 ) ~ (x) 2 = 0.07, 



(3) 



where (f(x)) — Y^k=i Pkf(xk)- Note that in the case of 
seismicity, Eq.© was found[23| to describe adequately 
the most probable value of II (a>). Furthermore, as shown 
in Ref.[23|], H(u>) for uj — > (or ki) can be considered as 
an order parameter for seismicity since its value changes 
abruptly when a main shock occurs and the statisti- 
cal properties of its fluctuations resemble those in other 
nonequilibrium systems (e.g., three-dimensional turbu- 
lent flow) as well as in equilibrium critical phenomena 
(e.g., two-dimensional Ising model). 

Apart from H(uj) or Ki, another useful quantity in nat- 
ural time is the entropy S, which is defined asJUEl] 
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S= (xlnx) - (x) ln(x). 
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This quantity depends on the sequential order of 
events [2^, l30ll and exhibits (3lj concavity, positivity and 
Lesche[33, [34j stability. The S value becomes equal to 
In 2/2 - 1/4 « 0.0966 for a "uniform" (u) distribution, 
as it was defined in Refs. [HI, m II H Hj, e.g. 
when all pk are equal or Qk are positive independent and 
identically distributed random variables of finite variance 
(In this case, k% and S are designated n u (— 1/12) and 
S u (= In 2/2 — 1/4), respectively). The same holds for 
the value of the entropy obtained [H, [3l| upon consid- 
ering the time reversal T, i.e., Tpk = PN-k+i, which 
is labelled by S-. The SES activities, when analyzed 
in natural time, have been found to exhibit infinitely 
ranged temporal correlations and -beyond Eq.Q- obey 
the conditions [2H I2H]: 

S,S^<S U . (4) 

The present paper is organized as follows: Section [TT1 
presents a new procedure for the study of the order pa- 
rameter of seismicity (/«i) in the area candidate to suffer 
a main shock during the period after the initiation of the 
SES activity until the mainshock occurrence. As exam- 
ples, we apply this procedure in Section IlIII to the most 
recent SES data. Finally, Section IIVI presents the con- 
clusions. 



II. NEW METHOD FOR SHORTENING THE 
TIME- WINDOW BY STUDYING THE ORDER 
PARAMETER OF SEISMICITY IN THE AREA 
CANDIDATE TO SUFFER A MAIN SHOCK 

Earthquakes exhibit complex correlations in space, 
time and magnitude (M), as shown by many studies, e.g . 

MJSJiE m, H3, m m M M M M M M m 

KM I4TJ . |41j . Moreover, it has been repeatedly proposed 
that the occurrence of earthquakes can be considered as 
a critical point (e.g. [42j], see also Ref.[43[ and references 
therein). Natural time reveals, as mentioned in Section I, 
that n(w) or k\ can be considered as an order parameter 
for seismicity [23|. To obtain II(cj) or Ki, however, it is 
necessary to decide the initiation time of seismicity anal- 
ysis. We decided to start the analysis immediately af- 
ter the SES initiation. This is based on our fundamental 
premise proposed long ago[l|,[ll|,[l2|] that both SES emis- 
sion and earthquake occurrence are critical phenomenon 
and, in a sense, the approach to "electrical" critical point 
shortly precedes "mechanical" critical point. 

A. Background 

Once the SES activity has been recorded, an estima- 
tion of the area to suffer a mainshock can be obtained on 
the basis of the so-called selectivity map [j. H3. H. |45| of 
the station at which the SES obsrvation was made. Thus, 
we have some area, hereafter labelled A, in which we 
count the small events (earthquakes) e.; that occur after 



the initiation of the SES activity. Each event et is charac- 
terized by its location x(ej), the conventional time of its 
occurrence t(a), and its magnitude M(e^) or the equiv- 
alent seismic moment Mo(ei). The index i = 1, 2, . . . , N 
increases by one each time a new earthquake with magni- 
tude M larger or equal to some threshold M t hres occurs 
within the area A (cf. Vi,i(ei+i) > t(ei)). Thus, a set 
of events A is formed until the mainshock occurs in A 
at i = N. To be more precise, a family of sets AM thr . es 
should have been formed before the mainshock occurs in 
A at i = N, where Au thres = {e> e A : M(e t ) > M thres } 
and the number of events in AM thrc3 is denoted by 
\-A-M thr e S |- The set AM thras becomes a (time) ordered 
set (i.e., the events are written in the sequence of their 
occurrence time) by selecting appropriately the indices 
j j for its elements e-,-, j — 1, 2, . . . , \AM thrc3 | so that 
Vj, t(ej + i) > t(ej). Since earthquakes do not occur ev- 
erywhere within the area A but in some specific loca- 
tions, we can also define R(AM thr ^ 3 ) as the minimal rect- 
angular (in latitude and longitude) region in which the 
epicenters of the events of AM thres are located. More- 
over, for a given ordered set AM thr£B the corresponding 
values of Ki(A Mthrss ), S(A Mthrs! ,) and 5_ (Awes ) can 
be obtained by analyzing in natural time its ordered el- 
ements ej (g AM thr<!S )- This is made by analyzing in 
natural time the pairs (Xj^Qj) = (i/l-^M t/1I . C3 1, Mo(ej)) 
where j = 1,2,..., |Ai th „J- 



B. The approach followed in our previous studies 

It has been repeatedly confirmed [HI, i i i 
l25l that, when an SES activity is observed, one 
can specify an area A by the selectivity map and (at 
least) one magnitude threshold Mthres which satisfy 
the conditions ([3|) and (|4]), i.e.,Ki(AM thras ) ~ 0.07 and 
S(^M tk „ s ), S- {AM thres ) < S u , a few days to a few days 
before the main shock. Thus, such a study enables, in 
principle, prediction of the main shock to be made within 
a few days to a few hours before its occurrence. 

The actual procedure was carried out as follows: For 
seismicity analysis, we set the natural time zero at the 
initiation of the SES, and then formed time series of 
seismic events in natural time each time when (small) 
EQs in A occurred, namely the number j increased. The 
normalized power spectrum in natural time for each of 
the time series was computed from the pairs (xk, Qk) — 
(k/j,M (ek)),k — l,...j, and compared with that of 
Eq.@ for (f> G [0, 0.5]. It was found that II(</>) approaches 
that of Eq. ^ a few hours to a few days before the main 
shock. We also calculated the evolution of the quanti- 
ties Ki, S and S- to ascertain Eq. (0} was also satisfied. 
The conditions for a true coincidence of observed time 
series with that of critical state were adopted as follows 
[H H llll, [Hill: First, the 'average' distance (D) 
between the H(4>) of the evolving seismicity and that of 
Eq.© should be smaller than 10~ 2 . This was a prac- 
tical criterion. Second, the final approach of K\ of the 
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evolving 11(0) to that of Eq.((2]), i.e., 0.07, must be by 
descending from above. This condition was found empir- 
ically. Third, both values S and 5_ should be smaller 
than S u at the coincidence. Finally and fourth, since the 
process concerned is self-similar (critical dynamics), the 
time of the occurrence of the (true) coincidence should 
not vary, in principle, upon changing (within reasonable 
limits) the the size of area and the magnitude threshold 

Mthres- 



C. A new approach suggested here upon using 

Vu thr( , s or the £[A M ] ensemble 

The basic idea behind the new approach suggested in 
this paper is the following: When area A reaches critical- 
ity, one expects in general that all its subareas have also 
reached criticality simultaneously. Each of these subar- 
eas corresponds to a certain value of Ki since the events 
that sequentially occur in that subarea (after the initi- 
ation of the SES activity) constitute a time-ordered set 
(which is a subset of -4M t; „. C3 ) , the natural time analysis 
of which results in a certain «i-value. Thus, we expect 
that the distribution of these Ki-values obtained from the 
analysis of all subareas should be centered around 0.070, 
according to Eq.Q. 

In principle, in order to investigate whether criticality 
has been approached immediately after the occurrence 
of a new event i within the predicted area A, we should 
construct all possible subsets that necessarily include the 
event i. Each of these subsets may be considered as a 
proper subset (VM thrss ) if and only if it includes all EQs 
that took place inside its corresponding rectangular sub- 
area, i.e., Ri^Pyithr^)- (This is, of course an approxi- 
mation -followed through out the present paper- because 
other geometries, e.g., circular, could be also considered.) 
In what follows, we will solely focus on such proper sub- 
sets of AM thrc3 ■ 

Let us now consider an example in which four earth- 
quakes shown in FigUJhave occurred in the area A, in a 
sequence indicated by the numbers 1,2,3 and 4. Figure [T] 
depicts the proper subsets of Ayi thrcB just after the occur- 
rence of each earthquake. In these proper subsets (which 
form the £[A M ] ensemble at each time instant), one 
has to compute the corresponding k% values and then 
construct their distribution Prob(Ki). The latter distri- 
bution depicts the state of the ensemble £[A M ]. For 
example, just after the occurrence of the second event a 
single proper subset can be defined, thus only ki[Ri(2)] 
is available. Just after the occurrence of the third event, 
three proper subsets of AM tllrss can be defined as shown 
in Fig. [TJ Recall that the necessary condition for having 
a proper subset at a given time instant is that it should 
include the last event (the third earthquake in this case) 
and for this reason the result corresponding to i?i(2) is 
not considered for the construction of the distribution 
Prob(Ki) [at this time instant]. By the same token, after 
the occurrence of the fourth event, seven proper subsets 



result, that can be visualized in Fig. [TJ Thus, we can 
now calculate m for each of these 7 subsets and after as- 
suming equipartition of probability among these subsets- 
we construct the Prob(Ki) versus K\ graph and then ex- 
amine whether it maximizes at k\ ks 0.070 (i.e., obeys 
Eq.©). 

In other words, in the new approach, the Ki-values of 
all these subareas that enclose earthquakes that occurred 
after the initiation of the SES until a given time instant 
are included together with the corresponding results of 
the largest area A, all of them treated on equal foot- 
ing. Then by performing an averaging procedure over all 
those proper subsets (which correspond, at a given time 
instant, to a large number of subareas of A), we expect 
that the average will also satisfy Eq. © . 

We shall demonstrate in the next section by using pre- 
cise experimental examples, which refer to the most re- 
cent SES activities recorded in Greece, that the following 
may also determine the time-window of the impending 
main shock: The average value {Ki('PM thr es)s[A M ] ob- 
tained when using the ensemble £[A M ] satisfies the 
condition © and then the mainshock occurs within a few 
days at the most. 



III. THE APPLICATION OF THE PROPOSED 
PROCEDURE TO THE MOST RECENT 
EXAMPLES 

In Fig[2J we depict four SES activities that have been 
recently recorded at the stations PAT (in central Greece) 
and PIR (western Greece): (a) on November 7, 2007, at 
PAT (b) on January 10, 2008, at PAT (c) on January 14, 
2008, at PIR and (d) on February 9, 2008, at PAT. These 
have been classified as SES activities by applying the 
usual criteria (he-, the conditions [3] and HJ) explained in 
detail in Ref . |46l. |47I. l48j . In addition, two long duration 
SES activities were recorded at PIR from January 21 to 
January 26, 2008 and from February 29 to March 2, 2008 
that will be described separately in Subsection E. We now 
apply the present procedure to all these cases: 



A. The case of the SES activity on Nov. 7, 2007 

The investigation of the seismicity subsequent to that 
SES activity was made in the area A:A r ff|£ , fo.o which 
was already indicated in Ref.[46j]. At 05:14 UT on 6th 
January, 2008, a major earthquake (M6.6) occurred with 
epicenter located at 37.1°N22.8°E, i.e., only around 50km 
to the south of the area A studied. What happened be- 
fore that EQ can be seen in Figj3] An inspection of 
this figure reveals that the Prob(fti) maximizes at around 
K\ = 0.073 and {ni{Vyi thTC3 ) s\A 1 = 0.070 with stan- 
dard deviation 0.008 upon the occurrence of a small event 
at 04:32 on January 4, 2008, i.e., almost two days before 
the main shock. 
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FIG. 1: The proper subsets and the corresponding rectangu- 
lar subareas immediately after the occurrence of the second 
earthquake "2" (upper panel), the third earthquake "3" (mid- 
dle panel) and the fourth earthquake "4" (bottom panel) . The 
location of each earthquake is shown by an open star. 



B. The case of the SES activity on Jan. 10, 2008 



The investigation of the seismicity was also made [47] 
in the area A-.N^^E^q q- The results are shown in FigEJ 
where we see that Prob(Ki) exhibits bimodal feature with 
a secondary peak at K\ 0.070 upon the occurrence of 
the small events at 12:21 UT and 13:26 UT on February 
3, 2008 (these two cases are shown with arrows). Ac- 
tually, at 20:25 UT and 22:15 UT on February 4, 2008 
two EQs with magnitudes 5.4 and 5.5 occurred with epi- 
centers around 38.1°N21.9°E lying at a small distance 
(ralOkm) from the measuring station PAT and inside the 
area studied. 



The case of the SES activity at PIR on Jan. 14, 
2008 



Since this SES activity was recorded at PIR, the study 
of the seismicity was made in the area -/V||.'q^|o.o m ~ 
dicated in Ref. [43] well in advance. The results of the 
computation are depicted in FigjSl which reveals that 




FIG. 2: The most recent SES activities recorded in Greece. 



the Prob(Ki) also exhibits bimodal feature, one mode of 
which maximizes at K\ f=a 0.07 upon the occurrence of a 
small event at 04:07 UT on February 12, 2008. Almost 
two days later, i.e., at 10:09 UT on February 14, 2008, 
a major EQ of magnitude 6.7 occurred at 36.5°N21.8°E. 
This EQ -according to USGS calatogue (which reported 
M w 6. 9)- is the strongest EQ that occurred in Greece dur- 
ing the last twenty years. In addition, a few hours later, 
i.e., at 12:08 UT, a M6.6 earthquake occurred at 36.2°N 
21.8°E. 
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D. The case of the SES activity at PAT on Feb. 9, 

2008 



The results until early in the morning on March 19, 
2008, of the study of the seismicity in the area [48 
-/VIy I-BIo.'o can be visualized in Fig. [6] which shows that 
Prob(Ki) has not exhibited a maximum at k\ w 0.07 yet. 
This investigation is still in progress. 



E. The case of the SES activity at PIR during the 
period February 29, 2008 to March 2, 2008 (Figf7](c)) 

This was a long duration SES activity (see Fig. Efc)) 
with polarity opposite to that of another long duration 
SES activity which was observed also at PIR from Jan- 
uary 21, 2008 to January 26, 2008 (FigJTJb)). Further, for 
the sake of comparison, in Fig[7{a), we present one more 
long duration SES activity on September 17, 2005 at PIR 
that was followed by the Mw6.7 EQ at 36.3°N23.2°E on 
January 8, 2006. The occurrence of the latter two major 
EQs (Jan.8, 2006 and Feb. 14, 2008), leads to the updated 
selectivity map of PIR as shown by the shaded area in 
FigJS] which lies along the Hellenic arc (marked with a 
thick solid line). In other words, the impending main 
shock for the SES activity of February 29, 2008 to March 
2, 2008 is expected to occur in this shaded area. 

A further study of the imminent seismicity is now in 
progress in order to clarify which region of the shaded 
area in FigJH will finally exhibit the expected feature 
,i.e., the maximization of Prob(Ki) at K\ «0.07 (see Ap- 
pendix) . 



IV. CONCLUSION 

Upon the recording of an SES activity, one can esti- 
mate an area A within which the impending mainshock is 
expected to occur. Following the subsequent seismicity, 
the probability density function of K\ is obtained, which 
maximizes at k\ w 0.070 a few days at the most before 
the occurrence of the mainshock. 



reported Mw=5.6), occurred at 35.0°N 25.3°E lying ap- 
proximately 150 km to the east of the southern part of 
the PIR selectivity map shown by the shaded area in 
FigH 

In view of the amplitude of the SES activity in Fig. 
[7f c) discussed in the main text (which is comparable to 
that of the SES activities depicted in Figs(7fa) and[7(b)), 
the study of the seismicity still continued to investigate if 
a more pronounced peak of Prob(Ki) at k,\ «0.070 would 
eventually occur in the near future. This actually hap- 
pened at 12:26 UT on May 8, 2008 and almost two days 
later, i.e., at 20:53 UT on May 10, 2008 a Ms(ATH)=5.6 
EQ occurred with an epicenter at 36.4°N22.3°E (see also 
below) . 



1. Note Added on May 19, 2008. 

In continuation of the study mentioned in the Ap- 
pendix, the ongoing investigation of the seismicity in the 
shaded area of Fig[5] reveals the following (cf. A calcula- 
tion with Mthres = 3.2 similar to the one in the previous 
cases, cannot be practically carried out in this case due 
to the large number of the events, i.e., more than half a 
thousand, involved in the calculation; recall that Mthres 
refer to the ML values reported by the Athens Observa- 
tory): 

(i) For Mthres — 3.8, a maximization of Prob(Ati) at 
Ki «0.070 was observed at 00:59 UT on May 16, 2008 
upon the occurrence of a ML=3.8 event at 36.1°N 21.8°E. 

(ii) For Mthres = 3.7, Prob(Ki) was clearly maximized 
at Ki «0.070 when a ML=3.8 event occurred at 36.6°N 
20.9°E at 20:27 UT on May 18, 2008. 

(iii) For Mthres = 3.6, Prob(Ki) exhibited a maximum 
at ki «0.070 at 10:28 UT on May 18, 2008 upon the 
occurrence of a ML=3.8 event at 35.9°N 23.3°E. 

The aforementioned calculations show that, at a first 
approximation and within reasonable time limits, mag- 
nitude threshold invariancc seems to hold. The extent to 
which, this behaviour will conform to the main conclu- 
sion of the present paper, it remains to be seen. 



2. Note Added on May 29, 2008. 



APPENDIX: EARTHQUAKES THAT FOLLOWED 
THE LONG DURATION SES ACTIVITY AT PIR 
FROM FEBRUARY 29, 2008 TO MARCH 2, 2008 

At 19:16 UT on March 25, 2008, the ongoing inves- 
tigation of the seismicity (for Mthres = 3.2) in the area 
Nfl.'o-Efo.'o (see Ref.jii]), after the SES activity of Febru- 
ary 29 to March 2, 2008 at PIR, showed a maximiza- 
tion of Prob(/ti) at K\ « 0.070 as shown in FigJ5] This 
was forwarded to various interested Institutes in Eu- 
rope, Japan and USA at 21:47 UT on March 26, 2008. 
Actually, almost one day later, i.e., at 00:16 UT on 
March 28, 2008, a Ms(ATH)=5.7 EQ (PDE of USGS 



At 15:11 UT on May 25, 2008 a Ms(ATH)=4.6 EQ 
occurred at 38.2°N 22.7°E (PDE of USGS reported 
mb=4.7). It is not yet clear whether this EQ is asso- 
ciated with the shaded area in FigJS] (in particular with 
its lobe that almost reaches PAT) or with the rectangular 
area N$$;gE$$. 

In addition, at 23:26 UT on May 27, 2008 a 
Ms(ATH)=5.1 EQ occurred with an epicenter around 
35.5°N 22.4°E, as expected by the Note added on May 
19, 2008. Upon the occurrence of this event, Prob(Ki) 
exhibits a pronounced maximum at n\ w 0.07 marked by 
an arrow in Fig llOl drawn for Mthres = 3.9. (An addi- 
tional arrow marks an earlier maximum on May 8, 2008 
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that preceded the aforementioned Ms(ATH)=5.6 EQ on 
May 10, 2008). Quite interestingly, this exhibits mag- 
nitude threshold invariance (a behavior that should be 
obeyed at the critical point) since a similar maximum 



at Ki = 0.07 appears simultaneously for M t hres — 4.0 
and M t h res = 4.1 as can be verified by an inspection of 
Figs. (fTT|) and (fl~2|) . respectively. 
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FIG. 7: The most recent SES activity recorded at PIR (c) along with the one that preceded the M m 6.9 EQ on February 14, 2008 (b) and the M w 6.7 EQ on January 
8, 2006 (a). 




FIG. 9: The arrow shows that the Prob(«i) of the seismicity (subsequent to the long duration SES activity recorded at PIR during February 29 to March 2, 2008) 
within the area iVff.'o -Efo.'o maximized at ki ~0.070 upon the occurrence of a small M3.2 event at 19:16 UT on March 25, 2008(i.e.,on the occurence of the 25th event 
after the SES). 




FIG. 10: Prob(Ki) versus Ki of the seismicity, for Mthres = 3.9, (subsequent to the long duration SES activity recorded at PIR during February 29 to March 2, 2008) 
within the shaded area shown in Fig(8] The two arrows mark the maxima at Ki = 0.07 that occurred on May 8, 2008 (i.e., on the occurence of the 29th event after the 
SES) and on May 27, 2008 (i.e., on the occurence of the 32nd event after the SES). The first maximum has been followed by the 5.6EQ on May 10, 2008, as described 
in the Appendix. 




FIG. 11: The same as Fie llOl but for Mthres = 4.0. The last histogram corresponds to the 5.1 event on May 27, 2008 and exhibits a maximum at Ki = 0.07. 



